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Abstract 

We present high resolution X-ray spectra of 14 cooling- 
flow X-ray clusters and groups obtained with the Reflec- 
tion Grating Spectrometers (RGS) on XMM-Newton. The 
spectra exhibit line emission from a number of Fe L charge 
states as well as O VIII, Mg XII, Ne X, Si XIV & XIII, N 
VII, and C VI. All spectra show a deficit of soft X-ray lines 
predicted from the isobaric multi-phase spectral model as 
compared with morphological mass deposition rates from 
spatially-resolved spectroscopy with the European Pho- 
ton Imaging Cameras (EPIC). We present some weak de- 
tections of plasma several times cooler than the ambient 
cluster temperatures. The results further suggest that ei- 
ther morphological mass deposition rates systematically 
overestimate the actual cooling rate or the emission mea- 
sure of cooling-flows has a much steeper distribution than 
that predicted by a simple isobaric multi-phase model. We 
briefly discuss some modifications to the cooling-flow pro- 
cess. 
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flow can be robustly predicted. If the blobs of gas cool in 
thermal isolation at constant pressure, and the dominant 
energy loss mechanism is via X-ray radiation, then the lu- 
minosity radiated per unit temperature interval must be 
proportional to the mass deposition rate, M: 



dL 
dT 
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where k is Boltzmann's constant, and \xm v is the mean 
molecular weight per particle. The only free parameter is 
M, which can be estimated from an X-ray image of the 
cluster core. The resulting spectrum can be calculated us- 
ing a collisional equilibrium spectral synthesis model with 
an assumed set of elemental abundances, and normalizing 
the contribution in each temperature interval as given in 
Equation (1). 

Data acquired by the Reflection Grating Spectrometer 
on XMM-Newton have enabled this robust spectral pre- 
diction to be quantitatively tested for the first time. Sur- 
prisingly, the observed spectra reveal a remarkable system- 
atic deficit of emission at low temperature, as compared 



to the multi-phase model (Peterson et al. 2001, Tamura 



ct al. 2001, Kaastra et al. 2001, Xu et al. 2002). Emission 



1. Introduction 

It has long been recognized that the cores of clusters of 
galaxies have sufficient X-ray luminosity to cool 10 to 



is observed at the expected levels for temperatures just 
below the ambient temperature of the cluster, T a (down 
to roughly \T a ), but not for significantly lower tempera- 
tures. This result has been confirmed with medium resolu- 
tion spatially-resolved spectroscopy using XMM-Newton 



1000 solar masses of keV plasma every year (e.g. Fabian yielded significantly smaller M's than expected ( 


David et 


fc Nuk-en 1977, Cowie & Binney 1977, Fabian 1994). The a l. 2001, 


Bbhringer et al. 2001, 


Molcndi & Pizzolato 2001, 


details of the cooling process are still debated, however, in [Schmidt, Allen, & Fabian 2001 


Ettori ct al. 2002). The in- 



most models parcels of cooling plasma collect at the center 
of the cluster, forming what is referred to as a cooling-flow. 
Direct evidence for cooler gas in the cores of clusters is in- 
deed wel l-established (|Canizares et al. 1979 , Canizarcs ct 



al. 1982, Mushotzky fe Szymkowiak 198? 



terpretation of this effect is still unclear, although a num- 
ber of possible physical mechanisms for suppressing the 
low temperature X-ray emission have been suggested (e.g. 
Peterson et al. 2001, Fabian et al. 2001). In this presen- 



If the gas were to cool homogeneously, the density pro- 
file of the cluster core would be relatively steep. That is 



inconsistent with imaging observations (lohnstone et al 



tation, we review the results of RGS observations for 14 
separate cooling-flow clusters, sampling a wide range in 
mass deposition rate. A companion analysis of the EPIC 
data for the identical sample of clusters is presented in 



1992), which has led to the conclusion that the cooling 
gas must condense locally in smaller clouds distributed 
over a large volume (tens of kpc), i.e. in a multi-phase 
medium. However, even ignoring the details of the result- 
ing spatial distribution, simple thermodynamic arguments 
show that the integrated X-ray spectrum of such a cooling 



Kaastra et al. (these proceedings) 



Diffuse X-ray Spectroscopy with the RGS 

We briefly comment on the spectral capabilities of the 
RGS, and the analysis of RGS data for extended sources. 



Proc. Symposium 'New Visions of the X-ray Universe in the XMM-Newton and Chandra Era' 26-30 November 2001, 
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Since, the RGS is a slitless spectrometer (den Herder ct 
al. 2000), its spectral resolution is degraded for extended 
sources. Nevertheless, in contrast to the transmission grat- 
ing instruments on Chandra, its dispersion is sufficiently 
high that it currently provides the highest spectral resolu- 
tion for the soft X-ray band for sources of arcminute size. 
The spectral resolution for an extended source larger than 
10 arcseconds is given by 



two line blends at 10.6 and 11.2 A roughly at the pre- 
dicted value of M. 

Previous, lower spectral resolution X-ray observations 
of cooling flow clusters had provided hints of this problem 
- a deviation from simple model predictions at the lowest 
energies. However, the discrepancies to model fits in those 
cases were generally interpreted as evidence for absorption 



AA sa ( b.lA (source size in arcminutes) / (spectral order)(2) 
At this resolution, one can still unambiguously resolve 
the most prominent resonance line contributions from all 
charge states of the most abundant elements. Of particular 
importance is the the Fe L complex, which can be resolved 
into its emission line contributions from Fe XXIV through 
Fe XVII (Li-like to Ne-like). The RGS resolution is several 
times higher than a CCD instrument in the Fe L region 
for a typical source size. 

The combined effective area of the two RGS instru- 
ments is 280 cm 2 . The field of view is 5 arcminutes in the 
cross-dispersion direction and 1 degree in the dispersion 
direction. This is sufficient to capture the entire cooling- 
flow region and produce a well-resolved spectrum with a 50 
ks exposure for a typical cluster. The instrument response 
for an oxtondod source is the convolution of the spatial die — 



by cool clouds embedded in the cooling flow ( White et al 
1991). For the RGS data, we can measure the intensities 



tributi bn of each omission lino and the off axis response for 

a given data selection region. To properly account for these 
effects, we use a novel Monte Carlo code for all spectral 



fitting analyses ( Peterson, Jernigan fc Kahn (in prepara- 
tion )|)~| A full spectral-spatial model of the source is used 
to set limits on the cooling-flow model. An example of the 
RGS data and corresponding EPIC-MOS data is shown in 
Figures 1 and 2. Figure 3 shows a Monte Carlo simulation 
for the same data set. 

The Soft X-ray Cooling-Flow Problem 

The failure of the isobaric multiphase cooling flow model 
to adequately describe the RGS data is illustrated for one 
of the most massive cooling flow sources, Abell 1835, in 
Figure 4. The predicted spectrum (at the resolution of 
the RGS for this source) is shown in the top panel. The 
model predicts detectable emission lines from the complete 
range of Fe L charge states, as well as from He-like and 
H-like charge states of N, O, Ne and Mg. The comparison 
of this model with the observed RGS spectrum of Abell 
1835 is shown in the second panel. As can be seen, the 
fit is reasonable for the higher ionization Fe L lines (Fe 
XXIV and Fe XXIII), but the model vastly overpredicts 
the emission line intensities for the lower charge states, 
especially Fe XVII. 

The final panel shows an empirical parameterization, 
where the isobaric model temperature distribution has 
been arbitrarily truncated at temperatures below kT = 
2.7 keV. This model has no clear physical motivation, but 
it does provide a good fit to the data. Note that it repro- 
duces the emission from Fe XXIV and Fe XXIII by the 



of the emission lines individually. It is clear from the raw 
data that absorption cannot explain this effect - the low 
energy lines are too weak, while the lower energy contin- 
uum is predicted correctly. Further, the RGS data make 
it clear that the failure of the model cannot be ascribed 
to the assumed elemental abundances. In particular, we 
see a deficit of low temperature emission lines in the Fe L 
spectrum alone. 

There is, however, clear detections of large quantities 
of plasma at temperatures above ^T a at roughly the pre- 
dicted quantity. This is very difficult to reconcile with any 
theoretical model. This problem is essentially distinct from 
the classic cooling-flow problem of the failure to detect 
the final products of cooling-flows, molecular cloud s and 
stars. In fact, recent detections of CO emission by Edge 



(2001) from massive cooling-flows now suggest that the 



end products of cooling may indeed be present. However, 
some unforeseen process must be suppressing the cooling 
radiation in the soft X-ray band. 

Other Clusters in the Sample 

Similar r esults for the full sample of 14 c lusters are pre- 
sented in Peterson et al. (in preparation). Generally, hot 



(several keV) clusters exhibit spectra similar to Abell 1835 
where no emission lines are detected below 3 keV. In- 
termediate temperature clusters (kT a « 3 keV) exhibit 
spectra similar to that shown in Figure 5 for the cluster 
2A0335+096. Prominent emission lines from Fe XXII-Fe 
XXIV are very apparent between 10 and 12 A. Weaker 
emission lines are also detected from the lower charge 
states, however they become increasingly more discrepant 
with the predictions as the temperature at which they are 
produced drops to lower and lower values. Since emission 
from essentially all Fe L ions is detected for these lower 
temperature clusters, the spectra clearly demonstrate the 
need for a continuous distribution of temperatures - but 
not the unique distribution predicted by the isobaric cool- 
ing flow models. 

A spectrum from a 1 keV cool group of galaxies, NGC 
533 is shown in Figure 6. Here again, there is evidence 
for significant cooling below ambient temperatures. The 
ratio of Fe XVII to Fe XVIII, however, is roughly 1, while 
the prediction for the standard isobaric model is closer 
to 3. Our sample includes 14 clusters, which differ in in- 
ferred cooling-flow mass deposition rate by three orders of 
magnitude. In every case, the isobaric cooling flow models 
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overpredicts the observed spectrum at the lowest temper- 
atures. It appears that the total soft X-ray luminosity is 
roughly consistent with the predicted morphological mass 
deposition rate, but that the emission measure distribu- 
tion is considerably steeper than this simple model would 



be liberated in Alfven waves or mildly relativistic particles 
as the magnetic field compresses. The very field lines that 
isolate the cooling blobs initially may reconnect and allow 
for modification of the cooling rate at each temperature. 



Finally, some heating models (e.g. by AGN outflows Pa- 



seem to require. 



bor & Binncy 1993, David ct al. 2001) could provide the 



Empirical Parameterization 

We are able to fit the observed spectra reasonably well, 
however, using a variety of empirical parametrizations. 
One that works especially well involves a simple modi- 
fication of Equation (1): 

dT 2 fiiiip T a 

Phis means, for example, that at temperatures of \ of the 
ambient temperature of the cluster, the X-ray luminosity 
is weaker than predicted by the model by a factor of 4 
or more. Phe integrated soft X-ray luminosity is therefore 
consistent with that needed for massive cooling-flows to 
within a factor of 2, but the results are clearly discrepant 
with the standard multi-phase model at the lowest tem- 
peratures. All clusters in the sample appear to follow this 
trend. 

Phe General Pheoretical Problem 

Phe cooling time for a constant pressure hot X-ray plasma 
is proportional to the temperature squared divided by the 
cooling function. Phis implies that it would take 6 times 
longer for Abell 1835 to cool from 8 keV to 3 keV than 
to cool from 3 to keV. It is extremely perplexing why 
there is no evidence for cooling after it has cooled 85% of 
the way. Phe lack of emission lines at lower temperatures 
could imply that either there is no gas cooling or it could 
also mean that the cooling rate at lower temperature is 
much faster than predicted. Phis could result from mixing 
processes which alter the rate of cooling at each tempera- 
ture or from additional coolants which carry energy away 
in addition to the X-radiation. 

Phe primary difficulty in finding a solution to this 
puzzle is that it is difficult to find a dynamical time- 
scale that would be so closely connected to the cooling 
time, which could vary by orders of magnitude depending 
on the local plasma conditions. Various ideas have been 
suggested, but none appears to naturally explain the ob- 



served phenomena. Mixing with cool clouds (Begelman & 



Fabian 1990) condensed from the cooling-flow could lib- 
erate some energy at other wavelengths, posssibly in the 



UV, as suggested by Ha observations (e.g. Heckman ct 
al. 1 981 |Crawford et al. 1999| ) and O VI flOegerle et al 
200lf)7or through PR dust emission (e.g. |Edge et al. 1999 . 



Allen et al. 2001). However, one would also naively expect 



necessary energy input to prevent the cooling, but they 
need to be finely tuned to work in just the right way for 
such a diverse sample of clusters. 

Phe solution to the soft X-ray problem could be com- 
plex and it could be difficult to distinguish between the al- 
ternatives. Future work connecting the X-ray observations 
to other wavelength observations, further empirical quan- 
tification of the temperature distributions in cooling-flows, 
and considerably more theoretical progress may result in 
a more complete understanding of the cooling process. 
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copious soft X-ray emission from the cloud interfaces with 
the hotter gas. Some fraction of the thermal energy could 
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Figure 1. An EPIC-MOS X-ray image for the cluster Abell S 1101. See Figure 2 for the corresponding RGS data. 
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Figure 2. The RGS data for Abell S 1101. Shown is the three images for the three possible projections of the data. 
The RGS data measures the dispersion coordinate, cross dispersion coordinate and CCD energy for each photon. The 
dispersion coordinate vs. CCD energy histogram shows the first and second order dispersed spectrum (curved lines) and 
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Isoboric Multiphase Cooling Flow Model 
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Abell 1835 and 2300 M Q yr 1 Cooling Flow + kT e =8 keV Ambient Component 
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Abell 1835 and Empirical Model with kT e =2.7 keV Cut Off 
4P — ' ' ' 1 ' ' ' ' 1 ' ' ' ' 1 ' ' ' ' 1 ' ' ' ' — =1 




Rest Wavelength (A) 

Figure 4- The three panels show the prediction for the isobaric cooling flow model, the comparison of the model with 
the spectrum of Abell 1835, and a model where the temperature distribution is cut-off below 2.7 keV. A full explanation 
is in the text. 
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Figure 5. The spectrum for the galaxy cluster 2A0335+096. The spectrum is typical of intermediate temperature (3 
keV) clusters. The cluster lies close to the galactic plane so there is significant absorption of the spectrum. Present in 
the spectrum is VIII (19 A), Mg XII (84 A), Si XIV & XIII (6.2 & 6.6 A), Fe XXIV -Fe XII (line blends between 
10.6 and 12.2 A, Fe XXI-Fe XIX (blends near 12.8 and 13.6 A), Fe XVIII (14.2 and 16.0 A), and Fe XVII (15 
and 17 A). The isobaric multi-phase model fails primarily because although low temperature ions such as Fe XVII are 
detected, they are much weaker than predicted. 
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Figure 6. The spectrum for the galaxy group, NGC 533. The spectrum is qualitatively similar to other low temperature 
(1 keV) systems. The line identifications are the same as the caption of figure 4- The isobaric multi-phase model fails 
in these systems primarily because Fe XVII is weaker than predicted relative to Fe XVIII. 



